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ABSTRACT 

We investigate the stellar structure of massive, quiescent galaxies at z ~ 2, based on Hubble Space 
Telescope /WFC3 imaging from the Early Release Science program. Our sample of 14 galaxies has 
stellar masses of M* > 10 10 8 M® and photometric redshifts of 1.5 < z < 2.5. In agreement with 
previous work, their half-light radii are < 2 kpc, much smaller than equally massive galaxies in the 
present-day universe. A significant subset of the sample appears highly flattened in projection, which 
implies, considering viewing angle statistics, that a significant fraction of the galaxies in our sample 
have pronounced disks. This is corroborated by two-dimensional surface brightness profile fits. We 
estimate that 65% ± 15% of the population of massive, quiescent z ~ 2 galaxies are disk-dominated. 
The median disk scale length is 1.5 kpc, substantially smaller than the disks of equally massive galaxies 
in the present-day universe. Our results provide strong observational evidence that the much-discussed 
ultra-dense high-redshift galaxies should generally be thought of as disk-like stellar systems with the 
majority of stars formed from gas that had time to settle into a disk. 



1. INTRODUCTION 

Structural evolution in the population of non- 
starforming (quiescent) galaxies has been observed and 
extensively discussed over the past few years. Such 
evolution has become an important ingredient in our 
description of evolutionary processes, and is used to 
constrain the theoretical framework for the formation 
of massive galaxies. The small sizes (i? e ff < 2 kpc) of 
massive (M ~ 10 11 M ) quiescent galaxies at z ~ 2 have 
received particular attention, both in terms of explaining 
their formation and their subse quent evolution (e.g., 
Truiillo et all l200i iDaddi et all l200l iPapovich etahl 
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The structural properties and surface brightness pro- 
files of these compact, massive z ~ 2 galaxies, beyond 
their small sizes, can help us understand the formation 
scenario for these remarkable objects and their subse- 
quent evolution: to produce such compact stellar sys- 
tems, highly dis sipative formation mechanisms have bee n 
proposed (e.g.. iHopkins et alll200l iWuvts et"alll2010f ). 
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Although it is not clear whether stars have time to 
settle into a disk in such a scenario, this may imply 
that these objects should be rotating. On the other 
hand, give n their high stellar masses , clustering prop- 
erties (e.g.. lOuadri et alll2007l: lHartlev et al.ll2010D . and 
number densities (e.g.. Ivan Dokkum et alll2010D . these 
objects must evolve into very massive, large, pressure- 
supported, bulge-dominated galaxies in the present-day 
universe. 

In this paper we explore the internal structures and 
surface brightness profiles of massive, quiescent galaxies 
at z ~ 2, utilizing Hubble Space Telescope (HST)/WFC3 
imaging, taken as part of the Early Release Science 
(ERS) program. We show that these galaxies are pre- 
dominantly disk-like and describe the differences with 
their likely present-day descendants. From this, a consis- 
tent narrative emerges in which highly dissipative events 
produce compact, disk-dominated stellar systems, which 
subsequently experience a series of merger events, simul- 
taneously explaining the growth in size and mass, and 
the lack of prominent disks in their descendants. 

2. SAMPLE SELECTION AND DATA DESCRIPTION 



FIREWORKS (|Wuvts et al.l 12008) provides a multi- 
wavelength catalog, from ground-based U band to 
Spitzer 24/im, and photometric redshifts. Stellar masses 
are derived for all sources with a minimum signal-to-noise 
ratio in the K band of 5, and photometric data points 
from at lea st four other bands. Fo llowing the procedure 
outlined bv lMarchesini et "all ([20091 ) . stellar-mass-to-light 
ratios are estimated f rom modeling the s p ectral energy 
distribution with the iBruzual fc Chariot) ([2003D stellar 
population synthesis model for solar metallicity and ex- 
ponentially de clining star formation rates. We adopt the 
iKroupal ([20011 ) stellar initial mass function, and the as- 
sumed cosmology is (12m, ^a, h) — (0.3, 0.7, 0.7). 

We select all galaxies with photometric redshifts 1.5 < 
z < 2.5 and stellar masses M > 10 10 8 M , at which the 
catalog is complete. Based on their rest-frame r — z and 
u — r colors we select those galaxies which are quiescent 
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Figure 1. Rest-frame u — r v s. r — z color for g alaxies drawn 
from the FIREWORKS catalog HWuvts et al.l 12003) . Dots repre- 
sent galaxies in the redshift range 1 < z < 2 and with stellar masses 
M > 10 10 M@, shown to demonstrate that galaxies separate into 
a star- forming sequence (below the dashed line), along which the 
exti nction increases towar d the top right, and a quiescent sequence 
(see [Williams ct al. 2009;, for a more clearly defined bi-modal dis- 
tribution). Large symbols represent galaxies with stellar masses 
M > 10 10 - 8 M Q in our redshift range of interest (1.5 < z < 2.5) 
and within the HST /WFC3 ERS mosaic. We select the 14 objects 
(indicated by filled circles) upward from the dashed line as our 
sample of massive, quiescent z ~ 2 galaxies. The objects indicated 
with open circles with crosses are considered star forming and not 
included in our sample. Two galaxies, indicated by open circles, 
that satisfy the color- criteria are excluded from the analysis (see 
the text). 



(see Figure 1): galaxies with blue u—r colors are excluded 
on the basis of the directly observed young stellar popu- 
lations; galaxies with red r — z colors are excluded based 
on the red color long- ward of the 4000A/Balmer break, 
indicative o f reddening by dust . Such a technique was 
first used bv lWuvts et al.l (I2007T) and, s u bsequ ently, for a 
much larger sample by [Williams et ail f|2009T ) . Whereas 
these authors use the U — V and V — J color combination, 
our u — r and r — z color combination is essentially equally 
effective (B. P. Holden et al. , 2011, in preparation). 

Sixteen of these galaxies are i n the Chandra 
Deep Field-South WFC3/ERS2 region (jWindhorst et all 
2010), for which an independent reduction of the 10 
pointings (5000s total) into an F160W mosaic was car- 
ried out by A. M. K oekemoer et al. (2011, in preparation) 
using Multidrizzle (Koeke moer et aDl2002h . 

Two of the 16 galaxies are excluded from further anal- 
ysis: one is a merger with strong tidal features; the other 
is located at a small projected distance from a very bright 
foreground galaxy, which precludes accurate photometry 
and analysis. 

We subject the images of the 14 remaining galaxies to 
the Lucy-Robertson deconvolution algorithm ("lucy" in 
IRAF), iterating 16 times, usin g a point sp read function 
(PSF) created with Tiny Tim (|Kristlll995f) . taking into 
account the dither pattern of the observations and the 



Multidrizzle data processing. The reconstruction brings 
out small-scale (^800 pc), high surface brightness fea- 
tures. 

We show the original and deconvolved images in Figure 
2. Structural parameters are inferred from fit ting single- 
comp onent Sersic models with GALFIT (P eng et alJ 
2002) to the original images, using the same PSF as de- 
scribed above. All 14 galaxies, which have a median stel- 
lar mass of 10 111 M®, have circularized half-light radii 
smaller than 2 kpc (the median is 1.2 kpc), in agreement 
with previous studies, a median Sersic index of 2.5, and 
a median axi s ratio of 0.67, alm ost identical to that of 
the sample of Ivan Dokkum et al.l 



We also determine, in the same manner, the structural 
parameters of the star-forming galaxies, that is, those 
below the dashed line in Figure 1. Their median Sersic 
index and axis ratio are n — 1.5 and b/a = 0.53, respec- 
tively. Their sizes vary with distance to the dashed line 
in Figure 1: those close to the line have sizes 1.5 — 2 kpc, 
comparable to the galaxies in our quiescent sample, while 
those far from the line are larger (4 — 5 kpc). This implies 
that our results do not depend on the exact criterion for 
quiescence. 

3. THE DISK-LIKE NATURE OF MASSIVE Z = 2 GALAXIES 

A significant fraction of the galaxies shown in Figure 
2 are strongly flattened in projection - the galaxies 2, 
7, 8, 11, and 12 have best-fitting axis ratios b/a < 0.5, 
indicative of a disk-like structure. Viewed more face-on, 
this would be difficult to recognize. Assuming approx- 
imately random viewing angles, this implies that a sig- 
nificant fraction - at face value a majority - of massive, 
quiescent z ~ 2 are disk-dominated. 

To follow up on the evidence based on the axis ratio 
distribution, we examine the disk-like nature of the in- 
dividual galaxies by performing two-component fits with 
GALFIT on the original images, using the same PSF as 
before. The positions of the two components are kept 
fixed at the position of the one-component fit, and the 
Sersic index is fixed at n = 1 for the disk- like component. 
All other parameters are left free. In Figure 3 we show 
4 examples of the two-component fits and point out the 
similarity between the two-component model fits and the 
deconvolved images shown in Figure 2. 

We explore the uncertainty in the two-component fits 
by artificially doubling the rms in the background by 
adding Gaussian noise, and by using two different PSF 
models: a PSF that is constructed from bright stars in 
the same mosaic, and a simpler TinyTim-based PSF. 
These alternative PSFs are known to be a less accurate 
representation of the true PSF than our preferred, de- 
fault PSF. 

For four out of 14 galaxies (nos. 2, 10, 11, 12) we re- 
produce the disk-dominated two-component fits regard- 
less of the adopted PSF model and the artificially in- 
creased noise level. For another three galaxies (nos. 4, 5, 
8) with disk-dominated two-component fits from our pre- 
ferred PSF and original noise level, the fitting parameters 
change with PSF choice and noise level to the extent that 
we consider the evidence for the disk-dominated nature 
of these objects as tentative. 

For these 4+3 galaxies either the % 2 value of the two- 
component fit is significantly lower than that of the one- 
component fit, or the one-component fit is disk- like itself 
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Figure 2. HST/WFC3 F160W imaging of our sample of massive, quiescent z ~ 2 galaxies (all with -Hfi60W < 23). The deconvolved 
images are shown directly below the original images. Photometric redshifts, stellar masses (in units of 10 11 Mq), Sersic indices, and axis 
ratios from one-component profile fits, as well as the scale of the images, are indicated. Ellipses indicate best-fitting axis ratios and sizes 
from GALFIT - the area of the ellipse corresponds to that of a circle with a radius that is twice the circularized half-light radius. All 
galaxies are clearly resolved and many are flattened in projection, indicative of a disk-like stellar structure. 



(with Sersic index n < 2). For the other 7 galaxies the 
fitting could not directly ascertain whether they are disk- 
dominated. 

To estimate the disk-dominated fraction among mas- 
sive, quiescent z ~ 2 galaxies, we combine the results 
from the quantitative two-component fits with the vi- 
sual impression reflected in Figure 2 (mostly flatness), 
which are independent of each other. All four galax- 
ies with robust disk-dominated two-component fits are 
flat in projection (with b/a < 0.6). We assign those 
objects weights of 0.9 (90% probability that these galax- 
ies are, in fact, disk-dominated). We assign weight 0.7 
to those 3 galaxies with more tentative quantitative ev- 
idence. These weights/probabilities should not be con- 
sidered precise and quantitative, but serve as a proxy 
for our level of confidence that these galaxies are disk- 
dominated. 

If we conservatively assume that none of the other 
seven galaxies are disk dominated (i.e., they have weight 
0), then we infer that 40% ± 15% of the population of 
massive, quiescent z ~ 2 galaxies is disk dominated. This 
number and its uncertainty include the weights as spec- 
ified above and the uncertainty due to the small sample 
size. 

However, some of the 7 non-classified galaxies, for ex- 
ample, 7 and 13, have small axis ratios. Therefore, it is 
not reasonable to interpret the lack of quantitative ev- 
idence for dominant disks as evidence for the absence 
of such disks. If we assign weight 0.5 to those objects, 
we infer that 65% ± 15% of the population of massive, 
quiescent z ~ 2 galaxies is disk-dominated. 

We have not attempted a more formal inversion of 



the intrinsic axis ratios as our sample is too small to 
de-project the observed axis ratio distribution without 
making overly restrictive assumptions about the under- 
lying intrinsic axis ratio distribution. Qualitatively, how- 
ever, the general flatness of the galaxies in our sample 
is consistent with a population with disks of finite in- 
trinsic thickness (~ 0.25). We stress that simply con- 
sidering the observation that 5 out of 14 galaxies have 
axis ratios b/a < 0.55 leads to a similar conclusion that 
~ 50% — 100% of the galaxies in the sample must be 
intrinsically thin, that is, disk-like: face-on counterparts 
of the observed edge-on disk-like galaxies must exist as 
well. 

We note that while very prolate systems can also lead 
to frequent small projected axis ratios, we consider this 
possibility unlikely: the absence of very prolate self- 
gravitating systems at lower redshifts suggests that na- 
ture does not produce such objects, and this should be 
independent of cosmic time. 

Based on examining the also available F125W and 
F098N imaging of our sample, we are confident that our 
inferences from the light profiles directly translate into 
information about the stellar mass distribution. The 
F125W and F160W images (roughly corresponding to 
B and V in the rest frame) are very similar in appear- 
ance, and all galaxies are much fainter in the bluer F098N 
passband (rest-frame U). We show color images for four 
disk-dominated galaxies in Figure 3; these illustrate that 
strong color gradients are absent, which implies that the 
disks are not strongly star forming^ In other words, 

8 None of these examples have 24/^m counterparts that signify 
star formation rates in excess of ~ 20MQyr~ 1 , which is less than 
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the light distribution must be similar to the stellar mass 
distribution. 

The major axis exponential scale length of the disks of 
the galaxies we consider as disk-dominated ranges from 
1 to 4 kpc, with a median of 1.5 kpc; this is substan- 
tially smaller than the scale lengths of disks in similarly 
massive galaxies in the present-day universe, which is 
- 4 kp c (log(iWkpc) = 0. 6 ± 0. 1 whe re 0.1 is the la 
scatter, iFathi et all (pOlOl ); iFathll (pOlOf )). The median 
scale length for our sample deviates by ~ 3tr, which im- 
plies that it is unlikely that most of these disk-like struc- 
tures will survive up until the present day. 

Several studies have before measured the size evolution 
of disk-l ike and bulge-like ga l axies based on their Sersic 
indices (|Truiillo et all I2006bt iBuitrago et ai]|2008t ). but 
only for stellar m ass selected samples t h at include star- 
formin g gala xies. iStockton et alJ (120061) , iMcGrath et alJ 



(2008), and IStockton et all (|2008t l reported the exis 
tence of passive, massive galaxies at z > 1.5 with disk- 
like surface brightness profiles, showing that a mix of 
morphological properties exists among passive, high- 
redshift galaxies. However, the popula tion o f com - 
pact galaxies reported by, f or example, iZirml (|2007f l 
and Ivan Dokkum et all (|2008l ) have not been identified 
with disk-dominated morphologies, but are generally 
thought of as early- type, bulge-like structures. Although 
Ivan Dokkum et all (|2008l ) mention the disk-like nature of 
some of the galaxies in their sample, they did not consider 
the implications for the overall population statistics. 

Our findings suggest that the majority of compact, 
massive, quiescent z ~ 2 galaxies are disk-dominated, 
typically hosting stellar disks with scale radii that are 
several times smaller than present-day disks. 

4. FORMATION AND EVOLUTION 

The scale lengths of the stellar disks in quiescent z ~ 2 
galaxies are, on average, ~ 3 times smaller than the scale 
lengths of similarly massive present-day stellar disks (see 
Section 3). This is consistent with the zeroth-order theo- 
retical expectation: the size and angular momentum of a 
z ~ 2 dark matter halo are on average ~ 3 times smaller 
than those of a similarly massive present-day dark mat- 
ter halo. The assumption here is that the correspondence 
between a disk and its halo in terms of re lative mass, size , 
and angular momentum does not evolve (|Mo et al.lll999l ). 
While these assumptions are obviously not strictly valid, 
it is still not unexpected that stellar disks are several 
times smaller at z ~ 2 than at the present. In addition, 
the global compactness of quiescent galaxies at z ~ 2 
is further boosted by the presence of compact, central 
components - these galaxies are not pure disks. 

Simulations of gas-rich mergers produce descendants 
that are broadly reminiscent of the observed galax- 
ies, most particularly in the sense that these simulated 
merger remnants have substantial rotatio n, and must 
therefore be disk- like ()Wuvts et al.l |2010() . Moreover, 
the simulations produce dense cores, and more extended 
outer parts. However, there are important quantitative 
differences. The merger remnants are entirely unlike 
exponential disks - they have very high Sersic indices, 
n > 4. This is the result of overly dominant central com- 

the average past star formation rate for these galaxies, qualifying 
them as quiescent. 



pact components. Gradual gas accretion to produce a 
disk, combined with the occasional (minor) merger event 
to produce a modest central component, may explain the 
described structural properties of quiescent z ~ 2 galax- 
ies. Generally speaking, their formation mechanism must 
be sufficiently gentle and slow as to allow gas to settle 
onto a disk before converting into stars. 

However, we note that the quiescent galaxies are not 
simply z ~ 2 star-forming galaxies of which the star for- 
mation rate has been reduced from 100 — 1000 Mq yr" 1 
(jForster Schreiber et all 12009( 1 to < 20 Mq yr" 1 .' Al- 
though massive star-forming galaxies at z ~ 2 are 
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20101) as well as in their dynamical propertie s (e.g., 
Forster Schreiber et al.ll2009t iGenzel et alJ 12010( 1. their 



disks are substantially larger (with scale lengths of ~ 
3 kpc) than the quiescent disks described here and 
only somewhat smaller than present-day disks. The 
explanation for this may be the age difference be- 
tween the quiescent and star-forming galaxies: the star- 
forming galaxies are undergoing a major, possibly first, 
episode of significant growth, while the quiescent galax- 
ies must have undergone such a p hase at an earlier 
epoch, implying smaller s izes (see also lFranx et alJ feOOS: 
Ivan der Wei et aT]|2009af) . The correlation between color 
and size for star-forming galaxies (Section 2) is consistent 
with this picture. 

Let us now consider the evolutionary path between 
the epoch of observation, z ~ 2, and the present day. 
In a hierarchical ACDM universe, the descendants of 
~ 10 11 M Q z ~ 2 galaxies are super-L* galaxies with 
masses > (2 — 3) x 10 11 Mq. This argument is bolstered 
by the observed clustering of massive z ~ 2 galaxies in 
comparison with the clustering properties of very mas- 
sive galaxies in the pres ent-day universe (jQuadri et alJ 
I2007t iHartlev et al.ll2010D . A simple, mo del- independent 
argument to the same effect is that the comoving number 
density of z ~ 2 galaxies with stellar masses ~ 10 11 Mq 
is the same as that of present-day galaxies with stellar 
masses - 3 x 10 11 M Q (jvan Dokkum et al.ll2010h . That 
is, these disk-like z ~ 2 galaxies are not progenitors of 
Milky Way type galaxies at the present day. 

Thus, we have compact, disk-dominated galaxies at 
z ~ 2, the present-day descendants of which are 2 — 3 
times more massive, are ~ 5 times larger in half-light 
radius, and almost never have prominent stellar disks 
(|van der Wei et al.l l2009bf ). Clearly, the structure of 
these galaxies has evolved dramatically over the past 
10 Gyr. Growth in mass through merging by the rel- 
atively modest factor of two or three compares wel l 
with observed merger ra tes (e.g.. iRobaina et ail 12010( 1 
and model expectations (|Hopkins et al.ll2010h . The ex- 
pected merger trees consist of a mix of frequent mi- 
nor mergers/accretion events and rare major mergers 
with low overall gas fractions (< 2 per galaxy since 
z = 2). The former are generally thought to provide 
the most efficient mechanism to explain size evolutio n 
(e.g., iBezanson et~aT1 l2009t Ivan der Wei et alJ l2009aT i. 
Indeed, the gradual buildup of the outer parts asso- 
ciated with such accreti on events has been observed 
van Dokkum et all 120101 ). To destroy a massive stel- 
lar disk, major merging is more efficient, although a se- 
quence of many minor accretion events can have the same 
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Figure 3. Top: F098N+F160W color composites for galaxies 11, 8, 2, and 5 from Figure 2, ordered by axis ratio. These four examples 
are chosen because of their flatness, with the exception of no. 5, which appears to have a compact bulge-like component surrounded by 
a more extended, disk-like component. There is no strong indications for color gradients, suggesting that the disk-components of these 
galaxies are not strongly star forming. Bottom: two-component model fits (without PSF smearing) for the same galaxies. The white and 
black ellipses indicate twice the size of the half-light ellipses for bulge-like and disk-like components, respectively. B/T is the ratio of the 
light in the model for the bulge-like component and the light of the models for the two components combined. "Rd" is the exponential scale 
length as measured along the major axis of the disk-like component, which we calculate by dividing the semimajor axis of the "half-light 
ellipse" by 1.6. "Rb,eff" is the circularized half-light radius of the bulge-like component. 



effect (jNaab et alJll999HBournaud et alJ f2007). Overall, 
a consistent narrative is emerging in which merging and 
accretion explain the growth in size and change in struc- 
ture of massive, passive galaxies over time. 

Further support for the direct link between the disk- 
like galaxies at z ~ 2 and the pressure-supported, 
massive elliptical galaxies in the present-day universe 
is provided by the comparable stellar densities of 
the z ~ 2 galaxies and the cores of present-day 
ellipticals dBezanson eFaLl [20091 iHopkins et ail [20091: 
Ivan Dokkum et alT 20101) Dynamical modeling of very 
massive nearby ellipticals has revealed that whereas the 
global rotation rate is small, the majority of the stars, 
even in the inner parts, are on disk-like orbits (see 
Ivan den Bosch et al.l 120081 , for an example). This may 
be the archaeological remnant of the disk-like nature of 
its progenitors. 



The authors thank the referee, Matthew Bershady, for 
helpful suggestions that helped improve the paper, and 
Kambiz Fathi for sharing his disk scale length measure- 
ments and comments on the manuscript. 
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